Polaronic slowing of fermionic impurities in lattice Bose-Fermi mixtures 
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We generalize the application of small polaron theory to ultracold gases of Ref. [2l[ to the case of 
Bose-Fermi mixtures, where both components are loaded into an optical lattice. In a suitable range of 
parameters, the mixture can be described within a Bogoliubov approach in the presence of fermionic 
(dynamic) impurities and an effective description in terms of polarons applies. In the dilute limit 
of the slow impurity regime, the hopping of fermionic particles is exponentially renormalized due to 
polaron formation, regardless of the sign of the Bose-Fermi interaction. This should lead to clear 
experimental signatures of polaronic effects, once the regime of interest is reached. The validity 
of our approach is analyzed in the light of currently available experiments. We provide results 
for the hopping renormalization factor for different values of temperature, density and Bose-Fermi 
interaction for three-dimensional 87 Rb — 40 K mixtures in optical lattice. 
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I. INTRODUCTION 

Polaron physics and, more in general, electron-phonon 
interactions are one of the most influential areas of mod- 
ern condensed matter physics 

and are believed to play a major role in the physics 
of High-T c superconductors [l|, Q and strongly correlated 
materials. 

Ultracold gases, on the other hand, allow for the inves- 
tigation of open issues in condensed matter using clean 
and highly tunable systems (see e.g. and Q)- In 

the context of polaron physics, the so-called spin-polaron, 
i.e. a single spin down impurity immersed in a Fermi 
sea of spin up particles has been realized as the extreme 
limit of imbalanced Fermi mixtures whenever N-tJNi 1 
and a remarkable agreement between theory (y, [7| and 
experiments [a, |9fl has been achieved. The original po- 
laron problem [lfj deals however with fermionic particles 
(electrons) interacting with lattice vibrations (phonons) , 
which obey the bosonic statistics, and is therefore some- 
how closer to the physics of Bose-Fermi mixtures. 

Bose-Fermi mixtures have been widely inve stig ated 
during the last few years both theoretically fn] - [l3l ] and 
experimentally [l44l6l |. The main focus however has 
been on the effect of the fermionic component of the mix- 
ture on the coherence properties of the bosonic conden- 
sate and on the supcrfluid-to-Mott insulator transition. 
In addition theoretical efforts were devoted to investi- 
gate the emergence of supersolid and other exotic phases 

MM- 

In a strongly imbalanced mixture of Nb bosons and Np 
spinless fermions with Nf/Nb <C 1, the dilute fermionic 
particles act as dynamic impurities in the bosonic con- 
densate. On the other hand, if one focuses on the 
fermionic component, the experimental setup closely re- 
sembles the polaronic problem in condensed matter, since 
fermionic atoms interact with phononic excitations of the 
condensate. Again, the main advantage of ultracold gases 
is that both the relative densities of the components and 
their mutual interactions can be tuned much more easily 



and to a larger extent than the corresponding condensed 
matter case. The extreme imbalanced limit allows for 
example to neglect the interactions between different po- 
larons and address the single polaron regime with relative 
simplicity. 

A remarkable achievement in this direction has been, 
e.g. , the recent experiment performed by the Bloch group 
fl6j . where a lattice Bose-Fermi mixture of 87 Rb — 40 K 
was studied, exploiting an interspecies Fano-Feshbach 
resonance to tune the Bose-Fermi scattering length clbf 
and varying the relative densities Np/Ns of the species. 
This allowed to study the effect of the interspecies inter- 
action and of the population imbalance between bosons 
and fermions on the transition from superfluid to Mott 
insulator in a very controlled way. 

A theoretical description of one-dimensional Bose- 
Fermi mixtures in terms of a Luttinger liquid of polarons 
has been proposed in [l7j . while the problem of polaron 
formation for a single impurity immersed in a homoge- 
neous condensate has been studied in (l8j and more re- 
cently in Other works addressed the emergence of 
polarons in the context of cold atoms in optical lattices 
j20L l2ll ] . They considered bosonic impurities loaded in 
an optical lattice and the whole system immersed in a 
large condensate of a different bosonic species. Only the 
impurities were affected by the lattice, allowing for an ar- 
bitrary slowing down of the impurities without perturb- 
ing the condensate (see Section 11111 for further details). 
This scenario could in principle be realized in experi- 
ments using a species-selective optical lattice. However 
this kind of setup, to our knowledge, has not yet been 
applied to Bose-Fermi mixtures, although several exper- 
imental schemes have been proposed [22[ and species- 
selective lattices have already been successfully applied 
to Bose-Bose mixtures 0, [24[ . In current experiments 
on Bose-Fermi mixtures in optical lattices, like e.g. in 
Ref. [l6|, both species are affected (though to a different 
extent) by the same optical lattice. In the latter case 
the tunneling properties of both species are intrinsically 
connected to each other, and the properties of the Bo- 
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goliubov modes of the condensate and their coupling to 
the fcrmionic particles are modified by the lattice. 

For these reasons, in this work we generalize the the- 
ory developed in Ref. [2l| to the case where both the 
fermions and the bosons move in the same optical lat- 
tice. Despite the presence of the lattice, we show that, in 
a suitable parameter regime, the bosonic condensate still 
sustains phonon-like excitations and the general frame- 
work developed for the homogeneous case in Ref. [2l[ 
still applies. We find that fcrmionic particles are expo- 
nentially slowed down by the interaction with the Bogoli- 
ubov modes of the condensate, due to polaron formation. 
We also discuss the relevance of our approach to current 
experiments on Bose-Fermi mixtures, analyzing the as- 
sumptions we made in order to derive our theory. For 
the specific case of a 87 Rb — 40 K mixture, we provide re- 
sults for the polaronic hopping renormalization of a sin- 
gle fermionic impurity in several experimental setups, i.e. 
for different values of the bosonic density, lattice depth, 
Bose-Fermi scattering length and temperature. This ef- 
fect can actually be measured by looking at the expan- 
sion of the fermionic component of the mixture in a lat- 
tice when the trapping potential is suddenly removed and 

n f /n b < i m. 

The layout of the paper is the following: in the next 
section we explain how under suitable conditions a Bose- 
Fermi mixture can be effectively described in terms of 
polarons and derive an expression for the fermionic hop- 
ping renormalization due to polaronic effects. In Section 
IIIII we analyze our assumptions within a generic experi- 
mental setup for Bose-Fermi mixtures in optical lattices. 
Results for the fermionic hopping renormalization in a 
87 Rb — 40 K mixture are provided in Section ITVl Conclu- 
sions are drawn in Section W\ 



II. THEORY 



investigation is assumed to be described by a single band 
Bose-Fermi Hubbard model, with the Hamiltonian 



H= Hb + Hp + Hbf where 

h b = -Jb @ih + h - c -) ~ i lB ^ B 

<i,j> 

£nf(nf-l) 

i 

-Jf £ (c\cj + h.c.) - [1fN f 



Hbf = U BF £ b\kc\ 



<i,j> 



(1) 

(2) 
(3) 
(4) 



The fermions do not interact (directly) with each other. 
However, as we will see later, they can still interact via 
boson-mediated interactions, a situation which is very 
similar to the BCS model for standard superconductiv- 
ity, where phonons mediate the attractive interactions 
between the electrons. We consider first the case Jf = 0, 
where the fermions act as a set of static impurities on 
the bosonic system and their position is specified by 
the (discrete) distribution function fi. We assume that 
the bosonic system in the presence of impurities can be 
treated within Bogoliubov theory, meaning that we con- 
sider the solution of the Gross-Pitaevskii (GP) equation 
in the presence of impurities and quantize the oscilla- 
tions around the classical deformed ground state. If the 
fermionic impurities were absent (fi = Vi or equiv- 
alently Ubf = 0), then for the unperturbed system of 
bosons we can write the following GP equation by ap- 
proximating the bosonic field operators with c-numbcrs 
% « V° and S] a ($•)• [US 



A. Gross-Pitaevskii theory for static impurities 

The derivation of this Section follows closely the one in 
Ref. 2l|, generalizing it to the case where single hyper- 
fine states of a bosonic and fcrmionic species are loaded 
together into an optical lattice generated by counter- 
propagating laser beams of wavelength A and frequency 
W£. For far off-resonant laser beams, the atoms expe- 
rience a potential Vb/f — Vb/f SiLi sir^^Xi/l) with 

V b/f = s B / F Er /F , where E? /F = 2 jf /F y 2 (h = 1) is 
the bosonic (fermionic) recoil energy, s B if denotes the 
dimensionlcss lattice depth for bosons and fermions in 
the respective recoil energy and I = A/2 is the lattice 
spacing, which we use as a unit length. We choose the 
fermionic recoil energy Ef7 as energy unit throughout 
the paper. 

The effect of the trapping potential is neglected and we 
postpone a thorough discussion about the correctness of 
our assumptions to the next section. The system under 



Jb ^ + U BB m 2 ^ = ^ 



(5) 



where the sum in the first term runs over the nearest 
neighbors of the lattice site i. In the case of a uniform 
system (in the lattice), this equation is trivially solved 
by ipi = y/no and hb = Ubb^ — zJb, where z is the 
coordination number of the lattice and no is the den- 
sity of particles per lattice site in the fully condensed 
state described by the classical GP theory (no quantum 
depletion of the condensate). In the presence of static 
fcrmionic impurities (Jp = 0), the previous result has to 
be modified because the condensate macroscopic wave- 
function is distorted by the impurities. If this distortion 
is sufficiently small then we can expand the classical field 
around the unperturbed solution, i.e. bi w ip® + Si and 
keep only the leading non-zero terms in the fluctuation 
5i . This approximation is valid if ^j- <C 1 and in this case 
the GP Hamiltonian has the form Hq p = Hq + Hs + Hu n 
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whore : 

Ho =-Jb E M)*Tp° + h.c] - /^ E \t 

<i,j> 
. UbB v-v i ,0i4 



I 2 



(6) 



H s =-Jb E + hx -} " E \ 5i 

<i,j> i 

+2U B bJ2\ S ^\ 2 



i 



(7) 

(8) 



where ft is the impurity distribution. The linear term 
in the fluctuation Si in Hg is identically zero because we 
choose ipf as the solution of the unperturbed GP equa- 
tion. 

By imposing the first derivative of this expression with 
respect to 6* to vanish and using the conditions ip® = 
^/no and —/is + Ubb^q = zJ, we obtain the following 
GP equation for the fluctuation field 6 

- Jb E + z JBS l + 2UBBnaSi + UBFy/no.fi = (9) 
This equation can be recast in the following form 

where we introduce the healing length of the condensate 



2J B 
Ubbtiq 



■I 



(11) 



Therefore the GP equation for the fluctuation field Si 
has the form of a discrete modified Hclmholtz equation 
where the impurity distribution acts as a source term. 
For a weakly interacting condensate the healing length £ 
is larger than the lattice spacing I and we can consider 
the continuum limit of Eq. (|10[) applying the same con- 
siderations discussed in [21|. The healing length £ fixes 
the typical scale for the variation in space of the fluctu- 
ation field Si due to the impurities. This means that the 
perturbation induced in the condensate by the impurities 
decays exponentially in space with the healing length £ 
and the condition ^ <C 1 (small perturbation of the 
condensate due to the impurities) implies (Ubb > 0) 



\Ubf\ 1 
Ubb n £, D 



IUbfUUbbtio) 
(2J B )« 



< 1 



(12) 



where D is the dimension of the system under consider- 
ation. Using the solution of Eq. (fTU|) . the GP Hamil- 
tonian provides the classical value of the ground state 
energy as a function of the impurity distribution fj, i.e. 
E = E d {fi). 



B. Bogoliubov corrections 

In the previous subsection we considered how the clas- 
sical condensate is distorted in the presence of static im- 
purities without including any quantum effects. We now 
consider the Bogoliubov excitations on top of the classi- 
cal theory, decomposing the bosonic quantum operators 
in a classical and a quantum part, i.e. i>, = ipi + Q\ where 
ipi = + Si is the solution of the GP equation described 
above. If we insert this expression in the Hamiltonian and 
retain only the terms up to the second order in the fluc- 
tuation fields, we obtain that all the linear terms in the 
fluctuation fields (classical and quantum) disappear since 
we have chosen the classical part as the solution of the 
GP equation in the presence of impurities. Therefore the 
Hamiltonian of the system has the form H = Hg+E cl (fi) 
where 



Hn 



Jb E + h -°~) - M-B E 

<i,j> i 

-2U BB J2 § 1 § ^\ 2 



Ubb 



E ( § M) 2 e\ + MM) (13) 



As evident from the expression above, the quantum part 
of the Hamiltonian is independent of the impurity dis- 
tribution and only depends on the unperturbed classical 
ground state through The quadratic Hamiltonian Hg 
can be diagonalizcd using a Bogoliubov transformation 
(see [28| for a general treatment) and we can re-express 
Hg in terms of the Bogoliubov modes of the condensate. 

This can be done introducing the following transforma- 
tion which express the original bosonic fluctuation oper- 
ators 9 in term of new bosonic operators /3 and j3\ i.e. 



E 

keFBZ 



A 



(14) 



where the sum runs over the quasimomenta k of the lat- 
tice within the first Brillouin zone (FBZ) and k = 
(the condensate) is excluded from the sum. To make 
Hg quadratic the coefficients u^i and Vk,i have the form 



1 =e lk - R -w k 



l 



vv = 



i e^t* (15) 



1 (16) 



where e£ = e k + zJb = 2J_b^^ 1 (1 — cos(fcj)) > 
is the (shifted) single particle spectrum in tight-binding 
approximation, hw^ = y/etjet + 2UBBn<o) is the energy 
of the Bogoliubov mode [2g, [27[ and N s is the number 
of lattice sites. The major difference in these expressions 
with respect to the continuum case treated in [2l| is that 
the FBZ provides a natural cut-off for the single-particle 
energy and therefore also for the energy of the Bogoliubov 
modes. In Fig. [T]wc sketched for comparison the shifted 
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FIG. 1. (a) Shifted single particle dispersion e£ and (b) Bo- 
goliubov spectrum fku^ in the k~ = plane for D — 3 and 
Jb = 0.029, Ubb = 0.11, Ubf = 0.065, n = 1. Energies 
are expressed in units of the fermionic recoil energy and 
momenta in units of l~ , where I is the lattice spacing. 



single-particle spectrum e k and the energy spectrum hujk 
of the Bogoliubov modes. 

Once expressed in terms of the Bogoliubov operators, 
the Hamiltonian of the system is diagonal and assumes 
the form 



H" 



+AE* + Hp Hp 



E 

keFBZ 



fiwk/3i/3] 



(17) 

where AE q is the quantum correction to the classical 
ground state energy due to the zero-point motion of Bo- 
goliubov modes. 

As already noticed in [2l| , since the Bogoliubov Hamil- 
tonian (|13|) does not depend on the impurities distribu- 
tion, this means that the Bogoliubov spectrum is unaf- 
fected by the position of the impurities and we have the 
same oscillation frequencies which we would have in the 
absence of the impurities. The equilibrium position of 
these condensate oscillations are however shifted by the 
presence of impurities. Since the Bogoliubov spectrum 
does not depend on the impurity positions, we can also 
switch the order of the steps in the preceding derivation 
and calculate the Bogoliubov theory around the unper- 
turbed ground state (no impurities), which is given by 
the assumption &, = ip® + 6i and keeping only terms up 
to second order in the fluctuation fields 0*. In this case we 
obtain the same expression as before for the Bogoliubov 
part while the classical part is not anymore the solution 



of GP equation in the presence of impurities. The Hamil- 
tonian operator now has the form H = Em> + Hg + Hu n , 
where 



H Hn = U BF ^^2f i [§}+§ i ] 



(18) 



and E^a is a c- number. Using the Bogoliubov transfor- 
mation (fT4"]) . we get in terms of the Bogoliubov modes 



H — Ejpo + Hg + Hu n 



(19) 
(20) 



Hg = J2 ^k/Sk 
keFBZ 

I 

inn =J2 E ^k[M iik /3 k + M* k /3 k ]/. (21) 



i keFBZ 



where 



M 



i,k 



UBF\fng 



k-Ri 



and 



= Ubf\ 



noe. 



(22) 



(23) 



In this case the new bosonic operators annihilate 
Bogoliubov excitations around the unperturbed ground 
state and therefore do not annihilate the real vacuum 
defined above which is distorted by the impurities. This 
leaves us with a generalized Holstein model with phonons 
coupled to the fermionic density [29|, [3(| • The main dif- 
ference with respect to the original Holstein model is that 
here we have a continuum of phonons instead of a single 
phononic mode. The high-energy phonon contribution 
in this expression is cut-off due to the FBZ, while in the 
the physical cut-off is provided by 
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continuum case 
the inverse of the typical localization length of the impu- 
rities in the Wannier states, which appears explicitly in 
the matrix elements M (see subsection III Dp . For static 
impurities this Hamiltonian can be diagonalized by intro- 
ducing a unitary Lang-Firsov [3l| transformation which 
shifts the equilibrium position of the condensate around 
the places where the impurities are localized 



U = exp 



j keFBZ 



U (24) 



This makes the Hamiltonian diagonal in the bosonic op- 
erators recovering Eq. (|17[) . 



C. Slowly moving impurities 

If Jf is not zero, i.e. if the impurities can move through 
the lattice, then the problem is not fully solved using 
(|24[) . However the Lang-Firsov transformation provides 
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physical insight on how to proceed. Indeed, introducing 
now the fermionic operators Ci and cj , we can repeat the 
steps above by simply replacing the density distribution 
fi with the density operator hi = c|c» everywhere. The 
Lang-Firsov transformation now acts simultaneously on 
fermionic and bosonic degrees of freedom 



U = cxp 



E E 

j keFBZ 



(25) 



Using the Bakcr-Hausdorff formula it is possible to show 
that tffffi = .' k V ; .l/,. k //,. UhjU^ = hj and 

UCjU^ = CjXj where the operator xj creates a coherent 
cloud of Bogoliubov modes around the position j, i.e. 



X] 



exp 



.keFBZ 



(26) 



The Lang-Firsov transformed Hamiltonian now has the 
form 

H L f= - Jf E (XiCi)HXjCj) - P-N F - ^ E 1 '•./"'" • 



0.020 




(Color online) Effective interaction potential Vij 



FIG. 2 

between the impurities for (Rj 



R i 



= and D = 3 



{Jb = 0.029, Ubb = 0.11, Ubf = 0.065, n = 1). The central 
peak for i = j is proportional to the polaron shift (V^j = 2E P ). 
Energies are expressed in units of and lengths in units of 
I. 



<hj> 



(27) 



keFBZ 



As already pointed out, the presence of Bogoliubov 
modes induces (offsite) interactions between the impu- 
rities. The interaction potential has the form 

Vi,j = 2 ^ hu k \M k \ 2 cos [k • (R< - Rj)] (28) 

keFBZ 

Therefore even though spin-polarized fermionic impu- 
rities do not interact directly with each other, their cou- 
pling to the Bogoliubov modes of the condensate creates 
an effective interaction between them. As shown in Fig. 
@, the interaction is always attractive (Vjj > 0) regard- 
less of the sign of the Bose-Fermi interaction and decays 
very fast with the distance between the impurities. More- 
over, during its motion, the impurity drags a Bogoliubov 
cloud and this affects its kinetic energy, as evident from 
the first term in Eq. f|2T[) . If the impurities are not mov- 
ing (Jf = 0), their energies are lowered by an amount of 
energy which represents the potential energy gain due to 
the interaction with the Bogoliubov cloud. This charac- 
teristic energy scale for static impurities is the polaron 
shift E p , where 



we expect that the impurity and its surrounding cloud 
will tunnel together like a composite object, i.e. they 
form a polaron. The expressions (|19I20I21|) and their 
Lang-Firsov transformed version (|27[) are in practice valid 
even beyond the condition (|30p (slow impurity regime), 
while the results we derive below assume that impurities 
are slow in the sense specified by Eq. (|50|) . In this sense 
this treatment is analogous to the the antiadiabatic limit 
(fast phonons) of the small polaron theory introduced by 
Holstein p9l . l30l | where the impurities move much slower 
than the time taken by the coherent cloud to rearrange 
itself. Here we would like to point out that in the expres- 
sion (|30[) , the bare fermionic hopping Jp is not compared 
with the bare bosonic hopping Jb but with E p which ex- 
plicitly depends also on the Bose-Bose interaction and 
the bosonic density. This is because in Bosc-Fcrmi mix- 
tures the excitations of the bosonic condensate and not 
the original bosonic particles play the role analogous to 
of phonons in the standard polaronic problem. 



D. Single fermionic impurity in the strong coupling 
small polaron regime 



E p = tujk\M k \ 2 = U 



2 

BF 



1 



E 



"oe k 



keFBZ 



s keFBZ y K/ 



(29) 



and ft = /Up + Ep in Eq. (|27|) . 

Whenever this energy scale is much larger than the 
hopping parameter of the impurities, i.e. 



c 



Jf 

En 



< 1 



(30) 



We consider now the case of a single fermionic impu- 
rity immersed in a much larger BEC. In real experiments 
this single polaron regime is realized whenever Nb 3> Np 
such that (n F )~ l l D » £//, where n F = N F /N S is the 
number of fermions per lattice site. This implies that 
the average interparticle distance is much larger than the 
healing length of the BEC, so that also interactions in- 
duced by the Bogoliubov modes can be neglected. For a 
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single impurity the Lang-Firsov Hamiltonian becomes 



H 1 . imp = J2 Mlft-* E (l i a i ) t (i j e j ))+£fl 



keFBZ 



<i,3> 



. (31) 

where Eq is a c-number. For Jp = the fermionic and 
bosonic part are completely disconnected and the im- 
purity can sit everywhere in the lattice with the same 
energy. If Jp is nonzero but ( = < 1 the polaron is 

the appropriate quasiparticlc and the hopping term can 
be treated clS cX small perturbation. We focus now on 
the regime of temperature fc^T -C E p , where incoher- 
ent phononic scattering is highly suppressed [2l|. The 
degeneracy of the Wannier states can be removed by in- 
troducing Bloch waves labeled by k' for the impurity and 
considering 



AE(k',{N k })-- 



(k', {N*}\-J P Y;(Xi%)KXj<h)\*, {^k}> 



<hJ> 



(32) 

where {N^} indicates the configuration of Bogoliubov 
modes. This matrix element can be calculated using 
standard techniques for phonons [Toj . If we assume ther- 
mally distributed phonons, we get that the bare hop- 
ping of the impurity Jp is exponentially renormalized to 
J'p = Jpe~ s and the renormalization factor S is given 
by 



S= 2^ |M„, k | 

keFBZ 



[l-cos(k-a)](27V k + i) 



(33) 



where 



Inserting the expression (|22p of 



e hu: k /k B T_ 1 ■ 

the matrix elements M, one obtains that S(T,Ubf) 
U BF f(T), where 



/CO 



1 



n e k 



^ (^k) 3 

keFBZ y K/ 



1 - cos(k • a)] (27V k (T) + 1) 



(34) 

In the practice, within this approach, the renormaliza- 
tion factor S is proportional to the square of the Bosc- 
Fermi interaction Ubf, while the factor / only depends 
on the condensate properties. As evident from the ex- 
pression above, S docs not depend on the sign of the 
Bosc-Fermi interaction. This results in a Gaussian de- 
pendence of the renormalized hopping on the Bose-Fermi 
interaction, i.e. J F = Jpe~ UsF f . The prefactor /, to- 
gether with its dependency on the temperature T and 
on other parameters like Ubb,Jb and no, can then be 
calculated independently and results will be presented in 
Section HV] for the 87 Rb - 40 K case. 

All the results presented in this Section are expressed 
in terms of the parameter no, the density of particles in 
the condensate. However a fixed value of no corresponds 
to different values of the bosonic density ub, whenever 
one changes the temperature T or the Hamiltonian pa- 
rameters Ubb, Jb- In order to compare data with exper- 
iments, the results have to be expressed in terms of the 



0.15 
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FIG. 3. (Color online) Sketch of the renormalization factor S 
as a function of Bose-Fermi coupling Ubf and temperature T 
(ks = 1) for fixed no. Energies are expressed in units of . 



bosonic density ns and no(T)| ns has to be calculated 
selfconsistcntly, as explained in the next section. This 
point is crucial for the understanding of the temperature 
dependence of the S factor. For example, for fixed values 
of no, S increases with temperature due to the increas- 
ing number of excited phonons, and therefore one would 
expect that the minimal slowing of the impurities occurs 
at T = 0, where only the zero point motion of the Bo- 
goliubov modes contributes. At the same time however, 
no decreases with T for fixed bosonic density ns and 
the overall temperature dependence is determined by the 
competition between thermal depletion of the condensate 
and thermal excitation of the Bogoliubov modes. The en- 
ergy spectrum of the Bogoliubov modes huj^ acquires a 
temperature dependence through no(T)\ nB which is miss- 
ing in the standard condensed matter case, where an in- 
creasing temperature only increases the phononic pop- 
ulation and therefore the S factor [3(|. This results in 
a mayor difference between the condensed matter and 
present case, and also suggests the existence of a richer 
temperature dependence of polaronic effects in the Bose- 
Fermi mixtures realization. Indeed, it is possible that 
different mixtures (or even the same mixture for different 
parameter range) show different slopes in S(T) or even 
a nonmonotonic behavior. A sketch of the dependence 
of the renormalization factor S on T and Ubf for fixed 
no is drawn in Fig. [3] This would be the relevant case 
whenever the depletion can be neglected in the range of 
parameters under investigation and no ~ tib- As shown 
in Section lTVl this is not the case for the 87 Rb— 40 K setup 
we considered. In any case it is worth noting that also at 
T = there is a sizable contribution to the S factor from 
all the Bogoliubov modes with quasimomenta within the 
FBZ. 

Whenever the condensate is unaffected by the lattice, 
e.g. in the case considered in 



211 ] , the Bogoliubov modes 



are labeled like plane waves and the sum runs over all 
the possible momenta q. In the lattice the FBZ provide 
a natural cut-off to the high-energy phonons contribu- 
tion, while the analogous role of the physical cutoff in 
the continuum case is played by an additional exponen- 
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tial decay of the matrix elements M for large momenta. 
Indeed the matrix elements M in the continuum case [2l| 
are given by 

M >«« \fiS fM) <35) 

where e q is the free-particle dispersion, Huiq is the energy 
of the Bogoliubov modes in the continuum pTj . 

/ 3 -(q) = J dr| X ,(r)| 2 cxp( iq .r) (36) 

and Xi( r ) 1S the Wannier wavefunction of the impurity 
localized at site j. Therefore M is proportional to the 
Fourier transform with respect to the momentum — q of 
the density profile of the impurity in a Wannier state. 
Using a Gaussian approximation of the Wannier wave- 
function, Eq. reduces to 

^.(^(expC-^Vexp^q-r) (37) 

where a is the width of the Gaussian wavefunction. 
Therefore the modulus of the matrix element M decays 
exponentially for q <C q c — l/er. In a deep lattice typi- 
cally a <C I and the cut-off is much larger than the Bril- 
louin zone. Therefore the continuum theory cannot be 
used to address quantitatively the lattice case, not even 
by introducing an effective mass m e to take the lattice 
in account. On the other hand a situation in which the 
bosonic condensate is unaffected by the lattice would be 
more favorable to realize the slow impurity regime de- 
scribed above, since the impurity can be slowed down 
arbitrarily without affecting the bosonic kinetic energy. 
In current experiments on Bose-Fcrmi mixtures in a lat- 
tice this is not the case as explained in the next section 
and the hopping parameters Jf and Jb are related to 
each other. 



III. LIMITS OF VALIDITY OF THE THEORY 

In this Section we analyze the assumptions we made 
in deriving our approach within a generic experimental 
setup for Bosc-Fcrmi mixtures in optical lattices. 

Since bosons and fermions are loaded into a single op- 
tical lattice of frequency wi , both species experience the 
same laser intensity I and V B , F = a B j f {ujj,)I , where 
a(wi) is the atomic dynamic polarizability at the laser 
frequency. We define the dimcnsionless parameter 

7 = — = -, s (38) 

which rules the ratio between the lattice depth experi- 
enced by bosons and fermions and therefore between their 
kinetic energies. For fixed atomic species, 7 can be varied 
by changing the wavelength of the optical lattice. 



In the presence of a interspecies Fano-Fcshbach res- 
onance, the Bose- Fermi scattering length asF can be 
tuned by a magnetic field B, while the Bose-Bose scatter- 
ing length a b can be considered in practice as constant in 
the same range of B. As already discussed previously, po- 
larized fermions do not interact directly with each other. 



In a real experiment, the atomic gas is confined in a 
trapping potential, which is not included in our approach. 
In the three-dimensional case and for sufficiently shallow 
trapping potentials, our results can still be considered 
locally in a LDA framework. The situation is very dif- 
ferent in D — 2. In this case the trap radically modifies 
the properties of the system, providing a cut-off to the 
long-wavelength Goldstonc modes, which would destroy 
the condensate at any finite temperature in the homo- 
geneous case. Therefore, even though formally our ap- 
proach to the homogeneous setup in the two-dimensional 
case is well defined at T — 0, the Bogoliubov treatment of 
the bosonic component immediately breaks down at any 
finite temperature in that case, and therefore we cannot 
directly apply our findings to a two-dimensional setup. 



In practice (quasi) two-dimensional systems are ob- 
tained by strongly increasing the optical lattice in one 
direction (e.g. in the z direction) such that the motion in 
this direction is frozen and only the zero point motion has 
to be considered [Hf. Under this condition it is possible 
to distinguish between two regimes. Whenever as <C cr z , 
where o~ z is the typical width of the onsite wavefunction 
in the direction of the tight confinement, then the scat- 
tering process is still essentially three-dimensional even 
if the motion is essentially two-dimensional. In this case 
the 3D scattering length can be safely used and the local 
interaction U increases with the confinement in the z di- 
rection since the overlap between the local wavefunctions 
is increased. For the case as < <r z a more refined treat- 
ment is required [33j]. In this case strong modifications 
of the interaction both in modulus and sign can occur in 
the system. The main feature of interest for the present 
paper is that the confinement can actually be used to 
further tune the interaction between the bosons in order 
to access different regimes in Bose-Fermi mixtures. Even 
though the existence of a real condensate has been pre- 
dicted for the case of trapped (quasi-) two-dimensional 
setups [H|, HH at low temperature, the application of the 
Bogoliubov approach necessarily needs the trapping po- 
tential to be explicitly included in the treatment. Since 
in this case a full numerical solution of the corresponding 
Bogoliubov theory is required, we postpone the analysis 
of this interesting case to the future despite its intrin- 
sic interest. On the other hand we expect the general 
conclusions of the paper to be still valid also in a 2D 
setup. For the reasons above, the approach is developed 
in the general case, when possible, but only results for 
the three-dimensional case are shown. 
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Single band Bose-Fermi Hubbard model 



Bogoliubov approach 



Our first assumption is that the Bose-Fermi mixture 
under investigation is described using the single-band 
Bose-Fermi Hubbard model ([T]). This requires that 
i) higher bands contributions, ii) non-local interaction 
terms, and iii) next-nearest neighbors hopping terms 
are negligible in the parameter range under investiga- 
tion. The first condition is particularly crucial for the 
bosonic component where the local density can take ar- 
bitrarily large values. For deep enough optical lattices, 
a Gaussian approximation can be used to estimate on- 
site parameters for our model. In this approximation 

Wj( x ) = rii=i W^fai) where W^{xi) is a Gaussian 
wavefunction in one dimension (xt = x,y,z for D = 3) 
localized around the site j of the optical lattice, i.e. 

W?{xi) = (7Ta 2 )- 1 ^e W [-(x i -R j ) 2 /(2a 2 )} (39) 



where <jb,f = \Jh/mB,F^B,F and hu>B,F 
2Er' F ^/SB,F- In the same approximation 

D/4 



Ubb = 



AirciB 
m B 



TT 2 S B 



Ubf 



2na 



D/2 



BF 



-1/2 . -1/2 



(40) 



(41) 



B T *F 

where we have set (E^ = H = I = 1) and m r = 
l/( m i? 1 + m F 1 )' anc ^ a <5-like pseudopotential has been 
used to model the interaction between particles. The 
Gaussian approximation provides a very poor estimate 
of the hopping parameter, which can be expressed in a 
simple way for deep enough lattices (large s) using the 
asymptotic solution of the Mathieu equation [36j j 



Jb/f/Et 



B/F 



4s 



3/4 
B/F 



cxp| 



iy/SB/Fl 



(42) 



In practice however this formula applies with reasonable 
accuracy only for s > 10, while for smaller s values the 
hopping term is overestimated and a direct numerical 
evaluation of the hopping parameter is required. Consis- 
tency with the model (JTJ) requires [TTJ] 



O-BFi &B cr_B,F <C I 



Ubb 



n B {n B - 1) < hio B (43) 



These conditions are reasonably well satisfied in prac- 
tice if the onsite density of bosons Ub is not too large. 
For increasing lattice depth s, the hopping parameter 
J decreases exponentially, while the interaction term is 
slightly increased because of the increasing onsite overlap 
of the Wannier orbitals. It is important to point out that 
since both species move in the same optical lattice, the 
following relation applies for large sf and fixed 7: 

J f m ^ 7 -(3/4) exp[(2 *^(^-i)] (44) 



Jb 



and therefore the ratio Jf / Jb is not constant for fixed 7 
but still depends on the lattice depth sf- 



Our approach is based on the possibility of describing 
the bosonic component of the mixture in the presence 
of static or slowly moving impurities within Bogoliubov 
approach. This requires in general that neither quantum 
nor thermal fluctuations are strong enough to substan- 
tially deplete the condensate, i.e. the condensate fraction 



4> = n /Nb < 1 



(45) 



needs to be close to 1. 

The parameter a introduced in Eq. (|12[) quantifies 
the effect of the impurities on the condensate wave- 
function, such that if a <C 1 we can expand the GP 
equation around the unperturbed solution in the ab- 
sence of impurities. It is worth mentioning that a is 
markedly dependent on the dimension of the system since 

= \U BF \(U BB n )i whne = \UbfI . e a ig indc _ 

(2.7b) 2 2j b 

pendent of the Bose-Bose interaction Ubb and the con- 
densate density no in D = 2. The condition a -C 1 
provides a constraint to the maximum value of the Bose- 
Fermi interaction where our theory can still safely ap- 
plied. Indeed we have 



I rjraax I 
\ u BF I 



I jjraax I 
\ u BF I 



2Jb 3 
(2J B )i 
(Ubbuo)'- 



in D 
in D 



(46) 



Strictly speaking, however, the condition a > 1 does 
not imply that the bosonic condensate cannot be de- 
scribed anymore within the Bogoliubov approach but 
only that the distortion of the condensate wavefunction 
due to the impurities is sizable and a full solution of the 
GP equation in the presence of impurities is required. 

In this sense we would expect that for \U\ > Ug l p x our 
theory still qualitatively applies, being however not any- 
more quantitatively accurate, if the condensate fraction 
4> of the mixture is close enough to 1. Whenever <f> is in- 
stead much smaller than 1, the Bogoliubov modes are not 
anymore the appropriate quasi-particles to describe the 
bosonic system and an alternative treatment is needed. 

As discussed in Section HH the Bogoliubov spectrum in 
our approach does not depend on the impurity distribu- 
tion and therefore the properties of the Bogoliubov modes 
can be estimated by applying the Bogoliubov approach 
for the pure system. The condensate density no is in gen- 
eral unknown and has to be calculated self-consistently 
within Bogoliubov theory for a given density n b and tem- 
perature T. This requires adding one more equation to 
our approach, i.e. the number equation of Bogoliubov 
theory in the condensed phase (2|| 



n B = n + 
1 / e£ + Ubbtio 



(47) 



E 

keFBZ 



N k (T) 



ej + Ubbtiq - huj^ 

2hjJ\r 
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which we solve numerically. As already discussed above, 
strictly in the homogeneous two-dimensional case, the 
Hohenberg-Mermin- Wagner theorem predicts that ther- 
mal fluctuations destroy the condensate for arbitrarily 
low-temperatures and Eq. (|4"T|) cannot be used for finite 
T in D = 2. 



£ parameter 

For a <C 1 and tp « 1 the condensate in the presence 
of static impurities can be safely described within the 
Bogoliubov approach presented in Section [TTJ Moreover 
within these approximations the Hamiltonian description 
given in Eq. ([27} applies also to the case of mobile im- 
purities. However, in order to obtain a simple expression 
for the renormalization factor S in the single impurity 
case, we had to assume that the fcrmionic hopping Jp is 
much smaller than the polaron shift E p (£ <C 1), where 
E p has to be calculated from the theory. According to 
the definition given in Eq. (|29l) , the polaron shift is given 
by E p = U BF g, where 



1 

N~ s 



V n o £ k 

keFBZ y KJ 



(48) 



and only depends on the properties of the bosonic compo- 
nent. Therefore since the polaron shift increases with the 
modulus of the Bose-Fermi interaction Ubf, the condi- 
tion f -C 1 limits the minimum value of the Bose-Fermi 
interaction for which the formalism can be applied to 

\U B F\~>\U^p\ = ^f 

Motivated by the large number of parameters present 
in the theory, we summarize the range of parameters 
where our approach can be applied by introducing the 
ratio 



R 



I U BF 
I jrmin 
U BF 



BF 
BF 



(49) 



such that for R > 1 there is a window of parameters 
where our approximations can be simultaneously satis- 
fied. Intuitively the condition ( < 1 requires that the 
fcrmionic impurities move much slower than the typi- 
cal time taken by the phononic cloud to rearrange itself. 
Therefore in general small values of the fcrmionic hop- 
ping Jp and 7 are more favorable to our approach, mean- 
ing that mixtures where the bosons move faster than the 
fcrmions would be in general a better choice to reach 
the regime under investigation, even though the energy 
scales involved also crucially depend on the interaction 
Ubb and on the density of the condensate. 



IV. HOPPING RENORMALIZATION IN A 
THREE-DIMENSIONAL 87 Rb 40 K MIXTURE 

To be more concrete in this Section, we refer to 
the most commonly studied Bose-Fermi mixture, i.e. a 



87 Rb— 40 K mixture loaded into an optical lattice in an ex- 
perimental setup similar to the one used in Ref. (l6| . For 
87 Rb — 40 K close to the interspecies Feshbach resonance, 
the bosonic scattering length is as ~ lOOao = 5.3nm, 
where eto is the Bohr radius. Increasing values of the 
Bose-Bose scattering length are generally unfavorable to 
our approach, since this increases both a and C (Eqs. (fT2]) 
and (|3"0)) ) and decreases the condensate fraction <j) (Eq. 
(|15J) ), if the other parameters stay unchanged. We found 
that smaller values of the Bose-Bose scattering length 
substantially enlarge the range of parameters where our 
approach is quantitatively valid. This suggest that dif- 
ferent mixtures with smaller ag, like the e Li — 23 Na mix- 
ture theoretically studied in [19j (ob ~ 53ao = 2,8nm 
for 23 Na |13|), could be even better candidates to real- 
ize the regime under investigation, if loaded into optical 
lattices. Since, however, experimental data about lattice 
Bose-Fermi mixtures involving those species are not yet 
available to the best of our knowledge, we decided in this 
work to focus on the 87 Rb — 40 K mixture and postpone 
the analysis of other mixtures to future publications. 

As explained in the previous section we only present 
results for the case D = 3, and we first focus on the T = 
case. Due to the large value of the Bose-Bose scattering 
length, we have to restrict our analysis to rather small 
values of the parameter 7, i.e. 7 = 1/3,1/2, where the 
bosons are substantially faster than the fcrmions [38| . 
Results are summarized in Fig. 01 

As evident in Figs. H£i andH}} for both setups the con- 
densate fraction (j> is relatively large and for 7 = 1/3 is 
always above 90%. Small values of the lattice depth are 
in general favorable to the consistence of Bogoliubov ap- 
proach since for a given as the bosonic system is less 
correlated. However sp cannot be reduced at will in or- 
der to stay in a parameter regime where the Hamiltonian 
in Eq. ([T]) applies. For 7=1/2 the quantum depletion 
of the condensate is much larger and only for the lowest 
lattice depth under investigation (sp = 10 — 12) the con- 
densate fraction is above 0.9. Considering the parameter 
R defined in Eq. (g§J) in Figs. Eh andgH, it is evident that 
only for 7 = 1/3 there is a sizable window of parameters 
where our approach is quantitatively valid, i.e. R > 1. 
We however realized that the strongest limitations arise 
from the condition a <C 1, which constrains the max- 
imum value of the Bose-Fermi interaction, rather than 
( < 1. As discussed in the previous section, whenever 
a <C 1 is violated but the Bogoliubov approach is still ex- 
pected to apply, we expect our results to be qualitatively 
valid and therefore show the results for 7 = 1/2 in the 
plot for comparison. Low bosonic densities and shallow 
lattices are favorable to the theory. Indeed, by reduc- 
ing the lattice depth sp, the range of densities where the 
theory applies is substantially increased. 



In Figs. Eh and Ef we plotted results for the parameter 
f(no(n.B), sf)\t=o defined in Eq. (j34]) . recalling that the 
renormalized hopping J F is related to / by the simple 
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FIG. 4. (Color online) (a-b) Condensate fraction cj>, (c,d) R parameter, and (e,f) / factor, defined in Eqs. (|45I49[) and (|34|) 
respectively, for two three-dimensional 87 Rb — 40 K setups with (a,c,e) 7 = 1/3 and (b,d,f) 7 = 1/2 as a function of the bosonic 
density ub and of the lattice depth sf at T = 0. / is expressed in units of (E^)~ 2 . The minimum value of the lattice depth 
sf is chosen such that sb > 5 in both cases. 



relation J F = Jpe~ u ^ F ^ . Both setups show a similar be- 
havior of the / parameter, though the numerical values 
are quite different for different setups and the data for 
7 = 1/2 are only shown for comparison since R < 1 in 
that case. For fixed (and small) values of the lattice depth 
sp, f increases quite rapidly with the density rig and 
then saturates. For larger values of sp instead, / shows 
a maximum at intermediate densities and then slightly 
decreases. However the decrease at large fig and sf and 



the non-monotonic behavior could result also from a loss 
of accuracy in our approach, since in that region of pa- 
rameters R < 1 in both setups. 

For finite temperature, we concentrate on a specific set 
of parameters, i.e. we choose 7 = 1/3, sf = 15 (sb = 5) 
and low-density ub = 0.1, which is the regime where 
our approach quantitatively applies in a larger parameter 
range. In a real experiment where the hopping renormal- 
ization is measured by observing the cloud expansion, the 
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FIG. 5. (Color online) Hopping renormalization J f /Jf for 
nj = 0.1,7 = 1/3 and sf = 15 as a function of the Bose- 
Fermi scattering length expressed in units of Bohr Radii ao 
and of the temperature T expressed in units of fermionic recoil 
energy Ef ■ 



initial configuration of the gas would be inhomogeneous 
due to the confining potential. This would mean that, 
once the trapping potential is removed, different parts of 
the fermionic cloud would expand in the lattice with a dif- 
ferent renormalized hopping J F due to local value of the 
bosonic density ub- We found that the renormalization 
factor S in general increases with increasing density Hb, 
which would imply that the edge of the cloud, where ub 
is smaller, will expand faster. The results shown in Fig. 
[5] for the hopping renormalization J T F / Jp as a function 
of the temperature T and of the Bose-Fermi scattering 
length asF arc representative of the experimental situ- 
ation at the edge of an expanding cloud and provide an 
upper bound for the renormalized hopping. The temper- 
ature range is chosen such that T < T max = E p (U^ x ), 
where E p is the polaron shift, and <j> > 0.9. We ex- 
pect our estimate to be quantitatively more accurate for 
40ao > \olbf\ > 60ao. 

As evident from Fig. [5j for fixed temperature the hop- 
ping renormalization as a function of clbf takes in our 
approach a Gaussian shape and the renormalized hopping 
decreases with the modulus of the Bose-Fermi scattering 
length, such that J F /Jp w 0.95 for |obf| ~ 50ao at the 
border of the bosonic cloud (ng =0.1). A much larger 
renormalization effect is expected for larger densities. 

The effect of the temperature is very small in the range 
of parameters investigated and it is hardy visible in Fig. 
[SJ We found that for fixed clbf the renormalization fac- 
tor S slightly decreases {Jp/ Jf increases) with increasing 
temperature, due to the dominant effect of the thermal 
depletion of the condensate. This trend is opposite to 
the one naively expected (see Fig. [3]) and also in con- 
trast with the condensed matter case, where S increases 
with T [30j . It is maybe worth mentioning that for higher 
densities (not shown) we found a nonmonotonic behavior 
in J F I Jp(T), due to the increasing relevance of the ther- 
mal population of Bogoliubov modes at larger T. Since, 
however, large densities are generally unfavorable to our 



approach, this effect could result as well from a loss of 
accuracy. 



V. CONCLUSIONS 

In this work we described the emergence of polaronic 
effects in Bose-Fermi mixtures in optical lattices. Our 
approach is closely related to Ref. [21| and is based on 
using the Bogoliubov approach to describe the bosonic 
component of the mixture, considering first static and 
then slowly moving fermionic impurities. The main dif- 
ference to the case addressed in Ref. [2lJ is that in our 
case both species are substantially affected by the same 
optical lattice, as in currently available experimental se- 
tups. We showed that the effect of the optical lattice 
on the bosonic condensate does not radically change the 
main conclusions for the homogeneous case [39| . However 
the range of experimental parameters where the approach 
applies is substantially modified, whenever fermions and 
bosons move through the same optical lattice, since their 
hopping parameters Jp and Jb are related to each other. 

For static impurities weakly coupled to the condensate 
(a <C 1), we have shown that an approximate treatment 
of the GP equation for the condensate wavefunction is 
possible and the Bogoliubov spectrum does not depend 
on the distribution of the fermionic impurities. However 
we expect the Bogoliubov treatment of the condensate in 
the presence of impurities to be still valid for larger a val- 
ues, provided the condensate fraction <j> is large. Conse- 
quently, our results are expected to be qualitatively valid 
even beyond the regime a <C 1. We postpone a more 
quantitative treatment of the the regime with larger a 
and (j> = 0(1), which would require a fully self-consistent 
GP + Bogoliubov approach, to a future work. 

We derived an analytical expression for the hopping 
renormalization of a single impurity in the regime C "C 1 
(slow impurity). This effect can be measured in ex- 
periments involving strongly imbalanced (Np/Ns <C 1) 
Bose-Fermi mixtures in optical lattice, by observing the 
expansion of the fermionic cloud in the lattice when the 
trapping potential is suddenly removed (25| . Within our 
approach, the fermionic hopping in this regime is ex- 
ponentially renormalized due to polaron formation, i.e. 
J' F = Jpe~ . In the relevant parameter range, the renor- 
malization factor S is found to be proportional to the 
square of the Bose-Fermi interaction. Therefore we ex- 
pect for J F a Gaussian dependence on the Bose-Fermi in- 
teraction Ubf (or equivalently agj?) and no dependence 
on the sign of Ubf {cibf)- This would lead to very strong 
experimental signatures of polaron physics once the con- 
sidered regime is reached. 

The temperature dependence of the renormalization 
factor S results from a competition between the thermal 
depletion of the condensate, which induces a temperature 
dependence on the Bogoliubov spectrum fiw^ through the 
condensate density no(T)\ nB , and the thermal population 
of phononic modes. This dependence of the phononic 



12 



spectrum on the temperature is missing in the standard 
condensed matter case, where the renormalization factor 
S always increases with T. 

In order to provide a better connection with experi- 
ments, we discussed the relevant parameter regime for a 
three-dimensional 87 Rb — 40 K mixture in an optical lat- 
tice. Due to the large value of the bosonic scattering 
length as, we have considered setups where Jp <C Jb 
and low bosonic densities. We found a sizable renormal- 
ization effect already for ng = 0.1. For a fixed value of 
the Bose- Fermi scattering length o,bf, the temperature 
dependence of the renormalized hopping is found to be 
opposite to the one naively expected. The renormaliza- 
tion factor S slightly decreases for increasing T, being 
dominated by the thermal depletion of the condensate 
rather than by the increasing phononic population. This 



represent a mayor difference with the standard condensed 
matter case. 



ACKNOWLEDGEMENTS 

We thank S. Will and I. Bloch for useful discussions 
and first pointing out the possibility of measuring pola- 
ronic effects in lattice Bose-Fcrmi mixtures. A. P. thanks 
S. Knoop for providing interesting details about optical 
lattices and the Physics Department of Lund University 
(Sweden) for its kind hospitality during the completion 
of this work. This work was supported by the Deutsche 
Forschungsgemeinschaft DFG via Sonderforschungsbere- 
ich SFB-TR/49 and Forschergruppe FOR 801. 



[1] A. Lanzara, P. V. Bogdanov, X. J. Zhou, S. A. Kellar, D. [19 
L. Feng, E. D. Lu, T. Yoshida, H. Eisaki, A. Fujimori, K. 
Kishio, J. -I. Shimoyama, T. Noda, S. Uchida, Z. Hussain, 
and Z. X. Shen, Nature 412, 510 (2001). [20 
J.T. Devreese and A.S. Alexandrov, Rep. Prog. Phys. 72, 
066501 (2009). [21 
R. Jordens, N. Strohmaier, K. Gunter, H. Moritz, and 
T.Esslinger, Nature 455, 204 (2008). [22 
U. Schneider, L. Hackermiiller, S. Will, Th. Best, I. 
Bloch, T. A. Costi, R. W. Helmes, D. Rasch, and A. [23 
Rosch, Science 322, 5907 (2008). 

G. -B. Jo, YR. Lee, J.H. Choi, C.A. Christensen, T.H. 
Kim, J.H. Thywissen, D.E. Pritchard, and W. Ketterle, [24 
Science 325, 1521 (2009). 

N. V. Prokof'ev and B. V. Svistunov, Phys. Rev. B 77, 
020408 (2008). [25 
C. Mora and F. Chevy, Phys. Rev. A 80 , 033607 (2009). [26 
A. Schirotzek, C. Wu, A. Sommer, and M. W. Zwierlein, 
Phys. Rev. Lett. 102, 230402 (2009). [27 
S. Nascimbene, N. Navon, K. Jiang, L. Tarruell, M. Te- 
ichmann, J. Mckeever, F. Chevy, and C. Salomon, Phys. [28 
Rev. Lett. 103, 170402 (2009). [29 
G. D. Mahan, "Many-particle Physics", Chapter 7 , Klu- [30 
ver Academic/ Plenum Publishers, third edition (2000). [31 
A. Albus, F. Illuminati, and J. Eisert, Phys. Rev. A 68, 
023606 (2003). [32 
I. Titvinidze, M. Snoek, and W. Hofstetter, Phys. Rev. 
Lett. 100, 100401 (2008). [33 
I. Titvinidze, M. Snoek, and W. Hofstetter, Phys. Rev. 
B 79, 144506 (2009). [34 
K. Giinther, T. Stoferle, H. Moritz, M. Kohl, and T. 
Esslinger, Phys. Rev. Lett. 96, 180402 (2006). [35 
S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, 
K. Sengstock, and K. Bongs, Phys. Rev. Lett. 96, 180403 [36 
(2006). 

Th. Best, S. Will, U. Schneider, L. Hackermiiller, D. van [37 
Oosten, I. Bloch, and D.-S. Liihmann, Phys. Rev. Lett. 
102, 030408 (2009). [38 
L. Mathey, D.-W. Wang, W. Hofstetter, M.D. Lukin, and 
Eugene Demler, Phys. Rev. Lett. 93, 120404 (2004). 
F. M. Ciucchetti and E. Timmermans, Phys. Rev. Lett. 
96, 210401 (2006). 



J. Tempere, W. Casteel, M. K. Oberthaler, S. Knoop, 
E. Timmermans, and J. T. Devreese, Phys. Rev. B 80, 
184504 (2009). 

M. Bruderer, A. Klein, SR. Clark, and D. Jaksch, Phys. 
Rev. A 76, 011605 (2007). 

M. Bruderer, A. Klein, SR. Clark, and D. Jaksch, New 
Journal of Physics 10, 033015 (2008). 
L. J. LeBlanc and J. H. Thywissen, Phys. Rev. A 75, 
053612 (2007). 

J. Catani, G. Barontini, G. Lamporesi, F. Rabatti, G. 

Thalhammer, F. Minardi, S. Stringari, and M. Inguscio, 

Phys. Rev. Lett. 103, 140401 (2009). 

G. Lamporesi, J. Catani, G. Barontini, Y. Nishida, M. 

Inguscio, and F. Minardi, Phys. Rev. Lett. 104, 153202 

(2010). 

S. Will and I. Bloch, "private communication". 

D. van Oosten, P. van der Straten, and H. T. C. Stoof, 

Phys. Rev. A 63, 053601 (2001). 

K. Burnett, M. Edwards, CW. Clark, and M. Shotter, 

J. Phys. B 35, 1671 (2002). 

A. L. Fetter, Annals of Physics, 70, 67 (1972). 

T. Holstein, Annals of Physics 8, 325 (1959). 

T. Holstein, Annals of Physics 8, 343 (1959). 

I. J. Lang and Y. A. Firsov, Zh, Eksp. Teor. Fiz. 43, 1843 

(1962). 

Z. Hadzibabic, P. Kriiger, M. Cheneau, S. P. Rath, and 
J. Dalibard, New J. Phys. 10, 045006 (2008). 
D. S. Petrov and G. V. Shlyapnikov, Phys. Rev. A 64, 
012706 (2001). 

D. S. Petrov, M. Holzmann, and G. V. Shlyapnikov, 
Phys. Rev. Lett. 84, 2551 (2000). 

C. Gies, B. van Zyl, S. A. Morgan, and D. A. W. Hutchin- 
son, Phys. Rev. A 69, 023616 (2004). 
W. Zwerger, J. Opt. B: Quantum Semiclass. Opt. 5, 9-16 
(2003). 

C. Samuelis, E. Tiesinga, T. Laue, M. Elbs, H. Knockcl, 
and. E. Tiemann, Phys. Rev. A 63, 012710 (2000). 
The condition 7 = 1/3, 1/2 can be realized by choosing 
A ?s 764, 762nm. These values of the lattice wavelength 
are however rather close to the atomic resonances, where 
a substantial heating effect from photon scattering is ex- 
pected. This could make harder the experimental obser- 



13 



vation on long timescales. 
] During the completion of this work, we become aware 
that similar expressions for the matrix elements M in 
Eq. (|22[) have been very recently derived independently 
in M. Bruderer, T. H. Johnson, S. R. Clark, D. Jaksch, A. 



Posazhennikova, and W. Belzig, arXiv: 1007.3828 for mix- 
tures of fermions and bosons in different optical lattices 
where tilting of the lattice for fermions was considered. 



